Static and Dynamic Phases for Vortex Matter with Attractive Interactions 
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Exotic vortex states with long range attraction and short range repulsion have recently been pro- 
posed to arise in superconducting hybrid structures and multi-band superconductors. Using large 
scale simulations we examine the static and dynamic properties of such vortex states interacting 
with random and periodic pinning. In the absence of pinning this system does not form patterns 
but instead completely phase separates. When pinning is present there is a transition from inhomo- 
geneous to homogeneous vortex configurations similar to a wetting phenomenon. Under an applied 
drive, a dynamical dewetting process can occur from a strongly pinned homogeneous state into 
pattern forming states. We show that a signature of the exotic vortex interactions under transport 
measurements is a robust double peak feature in the differential conductivity curves. 

PACS numbers: 74.25.Wx,74.25.Uv,89.75.Kd 



Vortices in superconductors interacting with random 
or periodic pinning provide a model system for study- 
ing the interplay between particle interactions that favor 
one type of ordering and substrate interactions that fa- 
vor a disordered state or some alternative symmetry ll|. 
In such systems a rich variety of liquid [1(, glassy |2j, 
and incommensurate static phases 0,0] arise, and under 
an external drive numerous dynamical phases [H 0-13] 
occur that produce distinct transport measurement fea- 
tures 0, fsMlOj] . Most vortex systems have strictly re- 
pulsive vortex-vortex interactions; however, modified in- 
teractions (MI) consisting of short range repulsion and 
long range attraction are proposed to occur in low-re su- 
perconductors Tl|, 13| or type-I/type-II hybrid supercon- 
ducting structures [14j . Recent imaging experiments in 
the multi-band superconductor MgB2 revealed inhomo- 
geneous, stripe-like vortex configurations 15, 1||, inter- 



preted as evidence that the multi-band nature of the 
system produces coexisting type-I and type-II vortices, 
giving long range attractive and short range repulsive 
vortex- vortex interactions ID, 15 , 17 - 27 1 . There has been 
some controversy regarding this interpretation due to the 
presence of quenched disorder or pinning, since even for 
vortices with purely repulsive interactions, strong pin- 
ning can produce highly inhomogeneous flux distribu- 
tions [m. 

Previous work has shown that a variety of periodic 
clump crystals and stripe structures can occur in systems 
with competing long range repulsive and short range at- 
tractive interactions or with two step repulsive interac- 
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tions [29H32IJ . It is not clear that such inhomogeneous 
structures would occur for interactions with long range 
attraction due to the low energy of a completely phase 
separated state. Numerical simulations of the proposed 
MI vortices found evidence of clump or labyrinth struc- 
tures 33, HJ; however, other studies of this system found 
complete phase separation |27j . In order to address the 
ground state ordering of the MI vortex system with and 



without quenched disorder, we perform large scale nu- 
merical simulations. We find that for long simulated an- 
nealing times and no pinning, pattern formation does not 
occur and the vortices completely phase separate into a 
single giant clump, while for shorter annealing times the 
system can become trapped in metastable states exhibit- 
ing some pattern formation. The addition of pinning 
induces a crossover from the phase separated state to a 
homogeneous, pattern- free vortex state. This behavior is 
very similar to the wetting-dewetting transition found for 
liquids of mutually attractive atoms on surfaces. When 
the atomic attraction dominates, the atoms clump or 
form a single dewetted drop, while when the surface at- 
traction dominates, the atoms spread out and wet the 
surface [35|. We find that under an applied drive, for 
strong pinning the driven states can form patterns such 
as stripes or labyrinths. We show that transport mea- 
surements provide a clear way to determine whether MI 
vortices exist in a sample, since the dynamic phase transi- 
tion into inhomogeneous clump or stripe states produces 
a robust double peak in the differential conductivity. Our 
results are general and can be applied to other systems of 
particles with competing long range attraction and short 
range repulsion in the presence of a substrate, such as col- 
loidal particles or vortices in Bose-Einstein condensates 

Simulation- We simulate a two-dimensional (2D) sys- 
tem of N v — 400 vortices and N p pinning sites with peri- 
odic boundary conditions in the x and y directions. The 
vortex dynamics are obtained by integrating the follow- 
ing equation of motion: riidRi/dt) = F™+Ff +F d +Ff , 
where r\ is the damping constant and Ri is the position of 
vortex i. We use the vortex- vortex interaction force pro- 
posed for the type-I/type-II hybrid materials and multi- 
band superconductors j2fl . |34] | , 
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Here -Ki(r) is the modified Bessel function, Rij = |Rj — 
Rj|, = (Rj — Rj)/Rij, and A and A/6 are the pen- 
etration depths in the two bands. In previous pin-free 
simulations of this model, the constants a and b were 
varied to change the relative strength of the attractive 
term [34]]. Long range attraction and short range repul- 
sion are obtained when a > b and b > 1.0. The coefficient 
a = Ki(r c )/Ki(br c ), and we set r c = 2.1A and b = 1.1. 
Since the interaction falls off rapidly at large distances we 
cut off the potential for > 8A. We model the pinning 
sites as attractive parabolic traps with maximum force 
F p and radius R p = 0.3A, as previously used in simula- 
tions of repulsive vortices 0, 0] ■ The vortex to pinning 
site ratio is measured in terms oiB^/B, where is the 
field at which there is one vortex per pinning site and 
B is the vortex density, held fixed throughout this work. 
The thermal force term has the following properties: 
(Ff(t)) = 0.0 and (F T '(f)if '(*')> = 2 V k B TS l3 S{t - t'). 
After performing simulated annealing we apply an ex- 
ternal drive = F^St and measure the resulting vortex 
velocity V x = N~ J2i " v i ' which would correspond 
to an experimentally measured current-voltage curve. 

We first study the model without pinning and con- 
duct simulated annealing studies by starting from a high 
temperature molten state and slowly cooling to T = 
during a time r a . For instantaneous annealing (small 
r n ) the vortices fall into a disordered assembly of clumps 
34], labyrinths, or voids. The structures coarsen as r Q is 
increased, and for long r a , the system completely phase 
separates into a single clump. Such complete phase sep- 
aration was observed in Landau-Ginsburg simulations of 
this same MI vortex model 27]. These results suggest 



that the stripes and inhomogeneous patterns imaged in 
MgB 2 [U, El, a weak pinning material, are not pro- 
duced by MI vortices. It is also possible that MI vortices 
are present in this material but that pinning arrests the 
phase separation. Even if pinning could stabilize pattern 
forming structures, it is also possible for inhomogeneous 
pinning to produce similar structures for vortices with 
only repulsive interactions, so the structures imaged in 
TBI 16 1 do not provide conclusive evidence of the nature 
of the vortex interaction potential. 

In Fig. 1 we show the vortex and pin positions from 
long r a simulations of samples with periodic and random 
pinning at B^/B = 0.4225. At F p = 0.3 for square pin- 
ning, Fig. 1(a) shows complete phase separation of the 
vortices into a single clump, with all pins empty except 
those directly under the clump. In Fig. 1(b) at F p = 0.9, 
each pin captures one vortex producing a homogeneous 
flux background, while the remaining vortices form a sin- 
gle clump. At F p = 1.5 in Fig. 1(c), each pin captures 
multiple vortices, and there are patches of unpinned in- 
terstitial vortices. For larger values of F p , all the vortices 
are trapped by pins. Similar behavior occurs for random 
pinning, as shown in Fig. l(d-f). The single clump state 
for weak pinning [Fig. 1(d)] is followed at higher pinning 



FIG. 1: The vortex positions (filled circles) and pinning 
sites (open circles) after annealing in the presence of periodic 
(a,b,c) or random (d,e,f) pinning for B^/B — 0.4225. (a) 
Dewetted (D) state at F p — 0.3 where most pinning sites 
are unoccupied, (b) Partially wetted (PW) state at F v = 
0.9. Each pin captures at least one vortex and the remaining 
vortices form a single clump, (c) Wetted (W) state at F p = 
1.5 where most vortices are trapped by pins, (d) D state at 
F p = 0.3. (e) PW state at F p = 0.9. (f) W state at F p = 1.5. 



strength by a phase with a coexisting clump and uniform 
flux background, shown in Fig. 1(e) for F p = 0.9. The 
clump decreases in size as F p increases, and eventually 
all the vortices become pinned as shown in Fig. 1(f) for 
F p = 1.5. We have investigated a range of values of F p 
and B^/B as well as other system parameters, and in 
general always find the same generic phases. We find no 
regimes where multiple clump, stripe, or labyrinth phases 
are stabilized, nor are there any phases that resemble the 
structures observed for systems with long range repul- 
sive and short range attractive interactions EB. In 
analogy to the spreading of liquid drops on surfaces [35| , 
we term the states illustrated in Fig. 1 dewetted (D), 
where all the vortices form a single clump; partially wet- 
ted (PW), where all the pins are occupied but there is 
still a vortex clump; and wetted (W), where the clump 
is absent and the majority of the vortices are trapped by 
pins. 

Figure 2 shows the locations of the D, PW, and W 
regimes in the F p versus B^j B phase diagram obtained 
for both periodic and random samples where we vary B^ 
by changing the pinning density. To identify the loca- 
tions of the phases we measure M n , the fraction of pins 
containing at least n vortices, as shown in the insets of 
Fig. 2. Both Mi and M 2 are low in the D phase; Mi is 
high and M2 is low in the PW phase; and Mi and M2 
are both high in the W state. We find that the system 
enters the W state at high F p and B^/B, while for weak 
F p or low B,p/B, the D state dominates. 

We now show that the most promising method for de- 
termining whether MI vortices are present in a sample 
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FIG. 2: F p vs B^/B phase diagram showing the D (O)j 
PW (+), and W (o) states for samples with fixed B and var- 
ied B$. (a) Periodic pinning, (b) Random pinning. Insets: 
M n , the fraction of pins containing at least n vortices, vs F p 
for periodic (upper inset) and random (lower inset) pinning. 
Squares: Mi; x's: M 2 . 



is by using an external drive. Studies of vortices with 
strictly repulsive interactions driven over random pin- 
ning find that a disordered pinned state can undergo a 
dynamical ordering transition in the moving state when 
the effective pinning is reduced The onset of the 

dynamic ordering is associated with a peak in the dV/dl 
curves in experiments @,H| or d(V x )/dFd curves in simu- 
lations [3,l9|, [lfj. For systems with MI vortices, strongly 
pinned samples form the uniformly pinned state illus- 
trated in Fig. l(c,f). A dynamical structural transition 
occurs under an applied drive due to the reduced effec- 
tiveness of the pinning in the moving state, which allows 
the attractive part of the vortex interactions to draw the 
vortices into a highly heterogeneous configuration. In 
Fig. 3 we show snapshots of the vortex positions in the 
driven state for random pinning samples with three dif- 
ferent pinning strengths at B^/B — 0.64. In the D state, 
two types of dynamics occur. When F p is very weak, the 
clump depins and remains intact; however, for higher F pi 
above depinning the clump sheds a trail of pinned vor- 
tices, as shown in Fig. 3(a) for F p = 0.3 and Fj = 0.05. 
The trailing edge of the clump becomes rarefied, leading 
to a decrease in the effective attraction between vortices 
that allows them to be torn away from the clump by 
the pins. As Fd increases, the trapped vortices depin 
and rejoin the clump, which retains its shape for higher 
drives. The reassembled clump for the F p = 0.3 sample 
is shown at Fd = 0.8 in Fig. 3(b). In Fig. 4(a) we plot 



FIG. 3: The vortex positions (filled circles) and pinning 
sites (open circles) in snapshots of the moving state for a 
sample with (a-f) random pinning and (g-i) periodic pinning 
at Bt/t/B — 0.64 under a drive F d ~k showing the formation of 
heterogeneous states at high drives, (a) F p = 0.3, Fd — 0.05 



(b) F p 
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0.3, F d = 0.8 (ii). (c) F p = 0.7, F d = 0.25 
= 0.7, F d = 3.0 (iv). (e) F p = 1.3, F d = 2.5 
1.3, F d = 3.8 (vi). (g) F p = 0.7, F d = 1.6 
vu). (h) F p = 1.3, F d = 0.96 (viii). (i) F p = 1.3, F d = 3.0 
ix). The labels i-vi [vii-ix] correspond to the drives marked 
in Fig. 4(a-c) [(e,f)]. 



d{V x ) I dFd and a v versus Fd for the random pinning sam- 
ple in Fig. 3(a). Here a v is the standard deviation of the 
vortex density calculated using a 2D grid with elements 
equal in length to the average pinning lattice constant a. 
Homogeneous vortex configurations give small values of 
a v , while a v is larger for heterogeneous configurations. 
The pinned clump state has a large value of a v . As Fd 
is increased, cr v initially drops when the moving clump 
sheds vortices, as in Fig. 3(a), and then increases again 
when the clump reforms for higher drives, as in Fig. 3(b). 
There is a weak two peak structure in d{V x )/dFd reflect- 
ing the two step depinning transition, with the second 
peak corresponding to the drive at which all of the vor- 
tices depin. Figure 3(c,d) shows a sample with F p = 0.7, 
which begins in a PW state. Above depinning, shown in 
Fig. 3(c) for Fd — 0.25, the clump disintegrates and the 
vortices flow in a meandering path, while at higher drives 
the vortices reorder into a much more heterogeneous state 
containing a clump and patches of filaments, as shown in 
Fig. 3(d) at Fd = 3.0. Figure 4(b) shows that <j v is high 
in the pinned heterogeneous state, passes through a lo- 
cal minimum near the first peak in d(V x )/dFd when the 
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FIG. 4: d{V x ) I dFd (left axes, dark blue curves) and a v (right 
axes, light green curves) vs Fd- Higher values of a v correspond 
to more heterogeneous vortex configurations, (a-c) A sample 
with random pinning at B^fB = 0.64 for (a) F v = 0.3, (b) 
F p = 0.7, and (c) F p — 1.3. The letters i-vi indicate the points 
at which the images in Fig. 3(a-f) were taken, (d-f) A sample 
with periodic pinning at B^/B — 0.64 for (d) F p = 0.3, (e) 
F p — 0.7, and (f) F p — 1.3. The letters vii-ix indicate the 
points at which the images in Fig. 3(g-i) were taken. 



clump breaks apart, and then increases again to a value 
higher than that at Fd = in the vicinity of the sec- 
ond peak in d(V x ) /dFd where the vortices form a moving 
clump and filament state. In Fig. 3(e,f) we illustrate the 
moving states at Fd — 2.5 and Fd = 3.8 for a sample 
with F p = 1.3 that has a W pinned state. The vortices 
form a moving stripe state near Fd — 2.75, and a portion 
of the stripes collapse back into clumps as Fd increases; 
however, even for high Fd, some stripes remain present. 
Figure 4(c) shows that the d(V x )/dFd curve has a pro- 
nounced double peak and that a v increases with increas- 
ing Fd as the vortex configuration becomes more homoge- 
neous. In general, for higher F p we find more stripe-like 
driven states, and the double peak feature in d(V x )/dFd 
becomes even more prominent. These results show that 
transport measures are the most useful approach for de- 
termining whether MI vortices are present in multiband 
or hybrid superconductors, since such vortices produce 
robust double peak features in the I-V curves. Com- 
plementary imaging experiments should show increasing 
heterogeneity of the vortex structure at higher drives. 
We note that the dynamical driven states found here are 
much more inhomogeneous than the states found for sys- 
tems with long range repulsive and short range attractive 



interactions [30j |. 

The dynamics for samples with periodic pinning show 
the same overall trends. In the D regime, the moving 
clump sheds a trail of vortices above depinning, and as 
Fd increases, all the vortices depin and the clump re- 
forms. The corresponding d(V x )/dFd in Fig. 4(d) has 
a double peak feature correlated with jumps in a v to 
higher values as the system becomes more heterogeneous 
at higher Fd- For the PW state, above depinning the 
clumps break apart and shed both pinned and intersti- 
tial vortices. As Fd increases, the clumps begin to re- 
form until all the vortices depin and form a coexisting 
clump and moving labyrinth state with vortices flowing 
along the pinning rows, as illustrated in Fig. 3(g). For 
Fd > 2.75, the labyrinth phase gradually collapses to an- 
other clump phase. Figure 4(e) shows the corresponding 
d{V x )/dFd and a v curves where the transition to a more 
heterogeneous state at higher Fd appears as an increase 
in a v . In the W state the initial depinning occurs via 
localized incommensurations flowing along the pinning 
rows in one-dimensional channels while the other vortices 
remain pinned, as shown in Fig. 3(h) for F p — 1.3 and 
Fd = 0.96. The vortices remain confined along the pin- 
ning rows above the second depinning transition, and as 
Fd increases, the moving rows transition into the moving 
labyrinth state shown in Fig. 3(i). Fig. 4(f) shows two 
pronounced peaks in d(V x )/dFd for this sample, while a v 
increases at higher drives. 

For systems with purely repulsive vortex interactions 
there is generally at most only weak hysteresis in the 
velocity-force curves. For MI vortices, we find hysteresis 
for the strong pinning cases with both random and pe- 
riodic pinning. The double peak in d(V x )/dFd is a very 
robust feature in these systems for strong pinning and 
persists for larger system sizes and various filling ratios 
B+/B. 

In summary, we have studied vortex matter with long 
range attraction and short range repulsion of the type 
proposed to occur in type-I/type-II hybrid structures and 
multiband superconductors. In the absence of pinning, 
this system does not form patterns or periodic arrays of 
stripes or bubbles, but instead completely phase sepa- 
rates. In the presence of random or periodic substrates 
we observe what we term a vortex wetting transition as 
a function of substrate strength or pinning density, with 
a transition from a heterogeneous to a homogeneous vor- 
tex state for increasing substrate strength. Our results 
show that transport measurements can be used to test 
for the existence of modified vortex interactions since, 
under an applied drive, the strongly pinned homogeneous 
states transition to highly heterogeneous states of moving 
stripes and clumps, producing robust double peak fea- 
tures in the differential conductivity. Our results are also 
relevant to the broader class of systems with long range 
attraction and short range repulsion on periodic and dis- 
ordered substrates, which can be realized in a number of 
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soft matter systems. 
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